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摘要 :植被 物候 学 作为 研究 植被 与 环境 条 件 相 互 作用 的 科学 ,在 全 球 气候 变化 的 大 背景 下 已 成 为 国际 热点 研究 领域 ,其 中 森林 
植被 在 调节 全 球 碳 平衡 维护 全 球 气候 稳定 的 过 程 中 有 着 至 关 重要 的 作用 。 随 着 遥感 技术 的 发 展 ,多 种 遥感 指数 被 应 用 到 森林 
植被 物候 研究 中 ,其 中 以 MODIS NDVI 和 EVI 应 用 最 为 广泛 ,而 叶绿素 灾 光 (SIF) 作 为 植被 光合 作用 的 “ 探 针 ”也 被 广泛 应 用 于 
森林 植被 物候 研究 中 。 为 了 探究 3 种 指数 在 森林 植被 物候 人 研究 中 的 差异 与 特性 ,本 文 以 长 白山 温带 红 松 阔 叶 林 通 量 观测 站 为 
研究 区 域 ,采用 模型 拟 合 结合 动态 冰 值 法 提取 2007—2013 条 林 物 候 特 征 参 数 ,并 使 用 通 量 数据 (总 初级 生产 力 CPP ) 进行 验证 。 
结果 表明 :NDVI 与 EVI SIF 相 比 ,表现 为 生长 季 开 始 时 间 与 结束 时 间 的 明显 提前 和 滞后 ,与 CPP 数据 偏差 较 大 , 且 夏 季 生 长 季 
峰 期 曲线 形态 过 宽 且 平坦 ,无 法 较 好 反映 生长 季 变 化 特征 ;EVI 相 较 于 NDVI 有 所 改善 ,整体 变化 趋势 与 SIF GPP 基本 吻合 ,但 
依然 存在 秋季 衰减 时 间 稍 迟 于 SIF 与 CPP 的 问题 ;SIF 昌 然 存在 夏季 又 降 现象 ,但 依然 与 GPP 数据 一 致 性 最 好 ,可 以 较 好 反映 
出 森林 植被 季节 变化 特征 。SIF 数据 与 植被 光合 作用 的 紧密 关联 使 其 在 植被 物候 研究 中 具有 优 于 植被 指数 的 准确 性 ,并 随 着 
遥感 平台 的 增加 和 反 演 方法 的 改善 ,将 会 在 多 尺度 、 多 类 型 的 植被 物候 监测 中 发 挥 更 加 重要 的 作用 。 
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Abstract: Vegetation phenology research has gained increasing attention because it is the best indicator of terrestrial 
ecosystem responses to climate change, and their consequences for ecosystem functioning. Forests play an important role in 
the terrestrial carbon cycle and maintain global climate stability. With the development of remote sensing technology , 


multiple remote sensing indices are applied to the study of forest vegetation phenology, among which Normalized Difference 


基金 项 目 : 国 家 重点 研发 计划 项 目 (2016YFB0501501,2017YFB0504000) ; 国家 自然 科学 基金 项 目 (41401110, 31400393) 
收 稿 日 期 :2017-08-21 ; 修订 日 期 ,2018-03- 20 
* 通讯 作者 Corresponding author. E-mail; zhoulei@ igsnrr.ac.cn 


http ://www.ecologica.cn 


10 期 刘 啸 添 ”等 :基于 多 种 遥感 植被 指数 .叶绿素 丈 光 与 CO, 通 量 数据 的 温 融 针 冰 混交 林 物 候 特 征 对 比分 析 3483 


Vegetation Index ( NDVI) and Enhanced Vegetation Index ( EVI ) derived from Moderate-Resolution Imaging 
Spectroradiometer( MODIS) are the most widely used. With the launch of Greenhouse Gases Observing Satellite( GOSAT) , 
Global Ozone Monitoring Experiment 2 ( GOME- 2) and Orbiting Carbon Observatory 2 ( OCO- 2) satellite, chlorophyll 
fluorescence( SIF) as a probe for photosynthesis of vegetation has been widely applied for studying the vegetation phenology 
in the global scope. In this study, we have analyzed the phenological characteristics of Pinus koraiensis and broad-leaved 
mixed forests in Changbaishan flux station from 2007 ~ 2013, using a double logistic function fitting and dynamic threshold 
method. Thereafter, the phenological characteristics, parameters, and time series curves of the three types of data were 
analyzed and compared. In addition, the validity of the results was confirmed by daily gross primary production( GPP ) from 
2007 to 2010. We found that time series of NDVI exhibits an earlier start of growth season ( SOS) date and a late end of 
growth season( EOS) date than that of EVI and SIF, and that the shape of the curve of growing season is too flat and broad 
to reflect the seasonal variations accurately owing to the saturation effect. The time series of EVI had. more pronounced 
seasonal characteristics, which was more consistent with GPP than that of NDVI, although the former showed a slightly later 
decline. SIF had the closest correlation with GPP and the best ability to track the seasonal cycle of photosynthesis and 
reflect the seasonal changes in forest growth. The close relationship of SIF data with photosynthesis indicated that SIF is 
more likely to play a better role in vegetation phenology monitoring than vegetation index..Moreover, the phenomenon of 
rapid decline in SIF around summer solstice, mentioned by many phenological studies, is consistent with the findings of the 
present study. Therefore, we have discussed the causes of this phenomenon from several aspects and provided a more 
reasonable explanation. Comparison of the three remote sensing indices in the present study suggested that SIF reflects 
seasonal variation in forest vegetation phenology better. With the increase in the number of remote sensing platforms and the 
improvement of inversion methods, SIF will play a more important role in multi-scale and multi-type vegetation phenology 


research. 


Key Words: vegetation phenology; vegetation Index; chlorophyll fluorescence; flux data; temperate mixed forest; dynamic 
threshold method 


植被 物候 学 是 研究 植物 (包括 农作物 ) 与 环境 条 件 (气候 水文. 土壤) 相互 作用 及 其 影响 机 制 的 学 科 … ， 
植被 的 生长 变化 特征 会 影响 生态 系统 的 结构 与 功能 ,如 叶片 面积 .光合 作用、 碳 循环 物种 构成 等 ” ,这 些 变 
化 反 过 来 又 会 影响 气候 系统 ,从 而 加 剧 气候 变化 。 在 众多 植被 类 型 中 ,和 森林 植被 占 地 约 1.4x10 万 hm , 占 全 
球面 积 的 9.4% ,储存 了 陆地 生态 系统 76% 一 98% 的 有 机 碳 , 在 调节 全 球 碳 平衡 .维护 全 球 气候 稳定 等 方面 具 
有 至 关 重 要 的 作用 三 ,因此 监测 森林 植被 物候 并 研究 其 周期 性 变化 规律 具有 重要 理论 意义 与 应 用 价值 。 

卫星 遥感 数据 能 够 在 多 种 尺度 上 为 各 类 生态 系统 提供 植被 物候 的 有 效 信息 ,因此 成 为 森林 物候 监测 
的 主要 手段 2 。 时 间 序 列 的 卫星 遥感 观测 数据 ,例如 MODIS 数据 被 广泛 应 用 于 提取 植被 物候 特征 的 研究 中 ， 
其 中 植被 指数 产品 能 够 在 较 大 时 空 尺 度 上 准确 反映 植被 绿 度 ,如 归 一 化 植被 指数 (NDVI) 与 增强 型 植被 指数 
CEVI) ,因此 广泛 应 用 于 植被 物候 的 研究 中 ”" 。 尽 管 MODIS 植被 指数 被 公认 为 是 质量 较 高 的 植被 指数 产品 
之 二 ,但 由 于 计算 公式 和 后 处 理 中 的 不 足 “ ,植被 指数 在 物候 监测 中 存在 着 较为 明显 的 缺陷 ,如 NDVI 2 
LLB SR 云雾 遮挡 影响 ,上 且 在 高 植被 覆盖 度 区 存在 易 饱 和 问题 EVE 虽然 改进 了 反 演 波段 的 计算 方 
法 使 其 对 植被 冠 层 结构 更 为 敏感 ,但 依然 存在 残留 气 溶 胶 的 噪声 影响 ,使 植被 物候 监测 结果 的 不 确定 性 大 大 
提高 。 

随 着 Greenhouse Gases Observing Satellite( GOSAT) ,Global Ozone Monitoring Experiment 2( GOME-2) ， 以 及 
Orbiting Carbon Observatory 2( OCO-2) 等 卫星 的 发 射 ,使 得 在 全 球 范 围 内 进行 SIF 数据 反 演 并 投入 到 植被 物候 
研究 中 成 为 可 能 。 叶 绿 素 效 光 作 为 光合 作用 的 副 产 物 ,不 局 限于 叶片 表 观 颜色 的 变化 观察 ,而 是 监测 植 
被 内 在 的 光合 作用 过 程 '" 。Joiner 等 ”采用 GOSAT SIF 与 MODIS EVI NDVI 数据 对 全 球 多 处 区 域 进行 植被 
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季节 变化 分 析 ,结果 表 明 SIF 与 植被 指数 具有 不 同 的 季节 性 以 及 相对 强度 变化 , 且 SIF 含有 独立 于 植被 指数 
的 信息 。 刘 新 杰 等 ' 中 利用 GOSAT 卫星 数据 进行 SIF 反 演 并 结合 MODIS EVI 对 我 国 西北 西南 .华北 地 区 进 
行 物候 研究 ,发现 SIF 时 间 序 列 曲线 较 EVI 季节 变化 更 加 明显 上 且 相 对 较 早 ,认为 SIF 数据 中 含有 EVI 数据 无 
法 表现 的 特殊 信息 。Walther 等 .小 采用 GOME-2 SIF 以 及 MODIS NDVI EVI 数据 ,对 中 高 纬 森林 进行 物候 对 
比 研究 ,发现 3 种 指数 的 季节 变化 曲线 存在 较 大 差异 。 由 此 可 见 ,NDVIEVI 与 SIF 数据 在 植被 物候 监测 上 差 
异 较 大 并 对 监测 结果 产生 显著 影响 。 

中 国 东 北 地 区 森林 覆盖 度 高 ,其 中 以 长 白山 地 区 森林 植被 种 类 最 为 丰富 且 季 节 特 征明 显 , 是 森林 物候 人 研 
究 的 理想 区 域 ” , 且 国内 鲜 有 利用 SIF 数据 对 该 地 区 森林 植被 进行 物候 监测 的 研究 ,因此 本 文采 用 目前 应 用 
最 广泛 的 植被 指数 MODIS NDVI,EVI 以 及 GOME-2 SIF 数据 ,选择 具有 中 高 纬 森林 代表 性 的 长 白山 温带 红 松 
益 叶 林 通 量 站 点 为 研究 区 域 ,对 3 种 指数 的 物候 特征 进行 对 比分 析 , 并 结合 通 量 站 观测 数据 估算 出 的 总 初级 
生产 力 (Gross Primary Productivity, GPP ) 数据 进行 验证 ,以 探索 3 种 指数 在 森林 物候 研究 中 的 具体 特性 。 


1 研究 站 点 概况 


长 白山 温带 红 松 阔 叶 林 通 量 观测 站 位 于 吉林 省 延边 朝鲜 族 自 治 州 安 图 县 王道 昌河 镇 。 该 站 地 处 长 白山 
目 然 保护 区 内 ,地 理 坐 标 42°24'9"N , 128°05'45"E ,海拔 高 度 738 m, 属 受 季 风 影 响 的 温 市 大 陆 性 气候 ,具有 显 
著 的 中 纬度 山地 气候 特征 ,春季 干旱 多 风 , 夏季 人 炎热 多 雨 ,冬季 干燥 寒冷 ,年 平均 气温 3.6% ,年 平均 降水 量 
713 mm, 主要 集中 在 6 一 8 H ,全 年 日 照 时 数 为 2271—2503 h ,无 霜 期 109—141 d, 

长 白山 的 植被 具有 典型 的 垂直 地 带 性 ,是 我 国 自然 生态 系统 保存 最 完整 的 地 区 之 一 ,是 拥有 大 量 物种 资 
源 的 生物 基因 库 。 通 量 塔下 垫 面 植被 为 阔 叶 红 松林 ,为 典型 的 地 带 性 植被 ,主要 建 群 树 种 有 红 松 ( Pinus 
koraiensis ) HA ( Tilia tuan) | 3 t BK ( Quercus mongolica) »7K HAA ( Fraxinus mandshurica Rupr)、 色 木 ( Acer 
mono Maxim) 等 ,平均 树 高 26 mi | 


2 数据 来 源 与 研究 方法 


2.1 数据 获取 与 预 处 理 
2.1.1 MODIS 植被 指数 

NDVI EVI 植被 指数 产品 采用 MOD13A2 V006 版 本 ,16 d 最 大 值 合成 ,每 年 23 期 数据 ,空间 分 辨 率 1 km, 
时 间 范 围 2007—2013 年 上 所 采用 的 MODIS 遥感 产品 均 来 自 LP DAAC 网 站 (https://lpdaacsve.cr. usgs. gov/ 
appeears/ ) 提供 的 AppEEARS( Application for Extracting and Exploring Analysis Ready Samples) 应 用 。 


TERT HH PR GEARS FE ai ts ES S 表 1 权 值 参照 表 
导致 数据 中 售 有 了 哄 声 影响 ,因此 利用 MOD13A2 中 与 Table 1 Weight reference 
一 期 数据 相对 应 的 “Data Pixel reliability” 数 据 集 ( 表 1) coe hoe 
KRH OK AS i m J et SC EXT NDVI, EVI 数 ; 人 
据 进 行 质量 控制 ; 虽 除 掉 数 值 标识 为 3 (ok Se ie BK 较 有 可 信 度 , 需 参 考 其 他 OA 文件 an 
云 遮挡 ) 质 量 差 的 数据 ,再 结合 质量 文件 对 数据 标识 为 2 数据 插 补 点 0.2 


1 的 数据 进行 处 理 , 史 除 掉 存 在 云 掩 等 噪声 影响 的 数 
据 , 最 后 利用 线性 插值 插 补 数据 并 赋予 权 值 。 

从 图 1 可 以 看 出 ,经 过 质量 控制 后 处 于 非 生 长 季 阶 段 (1 一 4 月 ,9 一 12 月 ) 的 数据 由 于 采集 条 件 较 差 , 数 
值 偏 低 被 吻 除 ,但 时 间 序 列 数 据 曲 线 波动 依然 较 大 ,存在 多 处 异常 值 ,虽然 质量 文件 中 未 显示 有 环境 条 件 的 干 
WAR ,但 仍 存在 传感器 精度 ,观测 角度 等 影响 ,因此 需要 进一步 去 除 噪声 影响 以 供 后 续 分 析 。 

2.1.2 ”叶绿素 殉 光 数据 
2007—2013 年 的 SIF 数据 采用 Kohler 等 利用 搭载 在 MetOp-A 卫星 上 的 GOME-2'7*! 测量 数据 中 波段 
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图 1 质量 控制 前 后 的 NDVI、EVI 时 间 序 列 曲线 
Fig.1 NDVI, EVI time series curve comparison before and after the quality contro 


NDVI:; 归 一 化 植被 指数 ,Normalized Difference Vegetation Index; EVI ; 4498 4! fH Ok 48 AL, Enhanced Vegetation Index 


范围 740 nm 的 数据 反 演 得 到 8 d 最 大 值 合成 ,每 年 46 WAH 2S TAIZ HE 0:5°x0.5° 的 数据 ,根据 长 白山 通 量 
站 点 坐标 提取 该 站 点 的 SIF 时 序数 据 并 利用 残 差 平方 和 进行 质量 控制 ,高手 2(m Wom sr nm ) 的 数据 被 
剔除 并 采用 线性 插值 搬 补 。 

SIF 数据 的 空间 分 辨 较 低 ,国内 外 对 于 物候 研究 中 SIF 数据 的 使 用 通常 将 其 他 数据 的 空间 分 辨 率 重 采样 
为 0.5°x0.5°% ”小 。 考虑 到 本 研究 是 基于 通 量 站 点 为 研究 区 域 ,数据 的 处 理 计算 均 以 站 点 为 中 心 ,区 域内 空间 
同 质 性 较 好 ,植被 类 型 一 致 ,参考 Joiner 等 物候 研究 成 果 , 该 研究 引入 SIF 相同 尺度 的 MPI-BGC 模型 
GPP 数据 作为 数据 误差 的 分 析 参 照 ,结果 表明 和 森林 站 点 的 SIE 数据 与 通 量 站 以 及 模型 CPP 值 的 一 致 性 均 较 
为 理想 ,空间 代表 性 差异 的 影响 存在 但 可 接受 , 故 在 本 文 申 为 了 尽 可 能 保证 遥感 数据 与 通 量 数据 的 空间 代表 
性 一 致 ,保留 了 遥感 数据 的 原始 空间 分 辨 率 。 
2.1.3 通 量 数据 

总 初级 生产 力 (GPP ) 是 植被 单位 时 间 内 通过 光合 作用 吸收 太阳 光 能 产生 有 机 物 的 总 量 ,又 称 总 第 一 性 
生产 力 1 ,能够 直接 且 准 确 地 反映 植被 生长 状态 。 

通 量 观测 以 较 高 的 时 间 频 率 连 续 测量 地 表 植 被 与 大 气 间 二 氧化 碳水 汽 以 及 能 量 的 交换 ,提供 了 植被 精 
确 的 代谢 活动 数据 ,高 度 连续 的 观测 数据 也 更 加 准确 地 反映 出 植被 生长 季 的 起 止 等 关键 物候 时 间 点 。 同 时 ， 
长 白山 温带 红 松 阔 叶 林 通 量 观测 站 数据 采集 源 区 基本 在 1 km 左右 ,与 季节 更 奉 和 大 气 条 件 稳定 与 否 相 关 ， 
总 体 上 与 MODIS 植被 指数 产品 空间 分 辩 率 一 致 ” 。 因 此 使 用 通 量 观测 直接 获取 的 净 生 态 系 统 碳 交换 量 
(NEE ) 并 估算 生态 系统 呼吸 (Re) 计算 出 的 GPP, 常 作为 “ 真 值 ” 检 验 壳 感 地 面 反 演 模型 的 结果 。 

由 于 通 量 数据 存在 观测 条 件 不 满足 相关 假设 与 传感器 标定 问题 ,涡流 的 随机 性 与 传感器 内 部 的 噪声 问题 
等 ,因此 需要 对 通 量 观测 数据 进行 坐标 轴 旋 转 、WPL 校正 和 储存 项 校正 以 消除 地 形 、 空 气 水 热 传 输 和 观测 高 
度 对 于 观测 数据 的 影响 PO! 。 由 于 夜间 大 气 层 结构 稳定 导致 夜间 观测 到 的 C0, 交换 通 量 存在 不 能 反映 真实 交 
换 通 量 的 可 能 性 ,因此 需要 对 夜间 汕 流 较 弱 情况 下 的 数据 进行 V* 殊 除 并 对 天 气 异 常 、 电 力 不 稳 定 以 及 虫 禽 
于 扰 等 造成 的 错误 数据 进行 剔除 。 

在 进行 数据 的 校正 和 插 补 后 ,将 每 日 每 半 小 时 间隔 的 连续 观测 数据 加 和 得 到 2003—2010 年 逐日 的 GPP 
时 间 序 列 数据 (图 2)。 

通 量 数据 在 本 文中 作为 验证 数据 ,被 当做 植被 物候 的 客观 “ 真 值 ” 曲线 ,由 于 时 间 分 辩 率 越 高 ,植被 的 物 
候 曲 线 会 保留 更 多 植被 生长 变化 的 细节 ,提供 更 符合 植被 真实 生长 状态 的 物候 信息 ,因此 并 未 对 GPP 进行 8 
d 的 最 大 值 合成 而 是 计算 成 每 日 的 通 量 数据 并 应 用 到 物候 曲线 分 析 中 。 
2.2 ”物候 特征 参数 提取 

本 研究 中 采用 模型 拟 合 结合 动态 冰 值 法 进行 物候 参数 提取 ,并 在 TIMESAT 软件 3.2 版 本 完成 ,选择 生长 
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总 初级 生产 力 GPP/(gC m? d’) 


2003 2004 2005 2006 2007 2008 2009 2010 
年 份 Year 


图 2 长 白山 站 2003 一 2010 年 GPP 季节 变动 曲线 图 
Fig.2 Seasonal variation of GPP in Changbai Mountain station for 2003—2010 


45 FF ua HY [A] (start of growth season, SOS) .生长 季 结 束 时 间 (end of growth season, EOS) ,生长 季 长 度 (length of 
growth season, LOS) ,季节 中 点 (time for the mid of the season) 4 个 参数 ,并 统计 时 序数 据 的 最 大 值 时 间 点 进行 
物候 分 析 。 

2.2.1 数据 寞 弟 值 处 理 

经 过 预 处 理 之 后 ,数据 仍 存 在 着 明显 的 异常 点 ,必然 会 影响 拟 合 后 的 曲线 形态 和 精确 度 ,因此 需要 利用 
TIMESAT 进行 异常 值 处 理 。 

本 文选 择 利用 STL( Seasonal and Trend decomposition using Loess ) 分 解法 结合 数据 权 值 文件 的 方式 对 3 种 
指数 进行 异常 值 处 理 。STL 分 解法 是 以 鲁 棒 局 部 加 权 回 归 帮 为 平滑 方法 的 时 间 序 列 分 解 方法 ,其 中 Loess 
(locally weighted scatterplot smoothing) 为 局 部 多 项 式 回 归 拟 合 * ,与 原始 权 值 文件 结合 使 用 可 更 好 地 考虑 到 
数据 特性 。 经 过 异常 值 处 理 后 3 种 指数 的 时 间 序 列 曲 线形 态 得 到 明显 改善 ,去 除 生长 季 与 非 生 长 季 的 异常 
值 。 由 于 SIF 是 每 年 46 期 数据 而 MODIS 植被 指数 为 每 年 23 期 数据 ,因此 曲线 形态 上 SIF 变化 更 加 细致 。 三 
种 指数 异 当 值 处 理 前 后 曲线 对 比如 图 3 所 示 , 红 色 虚 线 为 处 理 后 别 除 挥 的 异常 值 部 分 。 


归 一 化 植被 指数 NDVI 
Normalized difference 
vegetation index 
cS 
iN 


0 
2007 2008 2009 2010 2011 2012 2013 
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图 3 3 种 指数 时 间 序 列 曲线 处 理 前 后 对 照 图 
Fig.3 Three time series curve before and after the processing comparison 


SIF :叶绿素 荧光 ， chlorophyll fluorescence 
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2.2.2” 拟 合 模型 选择 与 物候 特征 参数 提取 

TIMESAT 3.2 中 共 提 供 了 3 种 拟 合 模型 ,考虑 到 模型 的 构建 原理 与 适用 性 ,针对 D-L 拟 合 ( Double logistic 
functions) 与 S-G 滤波 (Adaptive Savitzky—Golay filtering) 通过 式 1 计算 Pearson 相关 系数 ,并 绘制 3 种 数据 生 
长 季 数 据点 的 残 差 箱 线 图 比较 两 种 方法 保持 原 时 间 序 列 曲线 整体 形态 和 高 质量 点 的 能 


5 (index,, — index, ) (index,, — index, ) 


1 


(1) 


JÈ (index; - index, ) y > (index,, — index, ) i 

式 中 ,n 为 一 年 数据 总 期 数 , index, 为 拟 合 后 的 第 i 期 数据 , index, 为 拟 合 后 均值 , index, 为 拟 合 前 第 过 期 数 
据 , index, 为 拟 合 前 数据 均值 。 

根据 图 4 和 图 5 所 示 ,D-L 拟 合 法 在 拟 合 前 后 曲线 CO D-LAM S-G 滤 波 法 
相关 系数 上 明显 更 具 优势 ,只 在 SIF 数据 上 上 略 低 于 S-G 
滤波 ,但 也 达到 了 0.955 ,其 它 两 种 数据 分 别 为 0.977 和 
0.966 ,很 好 地 保留 了 原 时 间 序 列 曲线 的 整体 特征 。 在 
生长 季 高 质量 数据 点 的 偏差 分 布 上 ,D-L 模型 处 理 后 整 
体 偏差 较 小 且 分 布 更 为 均匀 集中 ,表明 在 较 好 地 保留 高 
质量 数据 点 的 同时 也 使 得 曲线 更 加 平滑 。 

经 过 对 两 种 拟 合 模型 对 比 后 ,选择 效果 更 为 理想 的 
D-L 模型 并 结合 动态 阔 值 法 ( 式 2) 提取 物候 特征 参数 。 


Pearson 相 关系 数 


Pearson correlation coefficient 


index, = (index,,,. — index,,,,) x 20% (2) 数据 类 型 Data type 
式 中 ,index, 代 表 该 指数 一 年 中 数据 最 大 值 ,indexii 代 图 4 3 种 指数 拟 合 前 后 曲线 相关 系数 


表 最 小 值 , 取 两 者 振幅 的 20% 作 为 生长 季 开 始 / 绪 束 的 | Fig.4 Correlation coefficients of three indices fitting curves 
国 值 ,index， 左 侧 曲线 高 于 index 时 即 为 进入 生长 季 :; 
右 侧 低 于 index, 则 为 生长 季 结 束 。 
指数 生长 季 中 点 时 间 根 据 曲 线 左 绒 和 右 端 分 别 增长 /减少 至 80% 幅度 的 时 间 点 取 均 值得 到 ( 式 3) ,并 统 
计 每 年 达到 生长 最 大 值 时 间 点 。 


tnig = average ( Li s0% 9! +,80% ) (3) 


式 中 ,ti 为 生长 季 中 起 时 间 ) ti s0% At, goo, 分 别 为 曲线 左 病 和 右 问 分 别 增长 /减少 至 80% 幅 度 的 时 间 操 。 
3 结果 


3.1 3 种 指数 时 间 序 列 曲线 整体 特征 

3 种 数据 的 7 年 时 间 序 列 曲线 呈现 明显 的 季节 性 变化 且 整 体 变化 一 致 (图 6) ,曲线 为 单 峰 形态 , 非 生 长 
季 阶 段 的 曲线 受 限 于 数据 质量 等 原因 数据 波动 明显 大 于 生长 季 , 且 冬季 曲线 数值 均 高 于 0。 

植被 指数 时 序 直 接 反 映 了 植被 绿 度 的 季节 性 变化 并 间接 表现 出 植物 光合 作用 的 强 弱 及 其 季节 和 年 际 差 
o NDVI 年 均值 在 0.51 一 0.57,7 年 数据 标准 差 为 0.24,EVI 年 均值 为 0.29 一 0.31,7 年 数据 标准 差 为 0.17 ， 
NDVI 相 较 于 EVI 整体 数值 偏 高 ,同期 数据 近 2 倍 于 EVI, NDVI 曲线 非 生长 季 存 在 较 大 波动 ,有 多 处 异常 峰 
值 ,而 EVI 在 相同 时 间 段 内 曲线 更 加 平滑 ,整体 分 布 更 加 集中 ,尤其 是 生长 季 阶 段 相 较 NDVI 曲线 波峰 较 宽 且 
平坦 的 形态 ,EVI 生长 季 波 峰 表现 为 “尖峰 ”形态 ,更 加 准确 反映 植被 生长 季 的 变化 特征 。 

NDVI 与 EVI 时 序曲 线 的 生长 季 形 态 大 致 对称 ,最 大 值 基本 在 生长 季 峰 期 的 中 间 位 置 ,而 SIF 时 序曲 每 一 
年 的 生长 季 最 大 值 时 间 要 早 于 生长 季 峰 期 的 时 间 中 点 , 较 早 达到 生长 峰值 ,之 后 迅速 衰减 进入 衰老 期 ,这 也 是 
SIF 数据 与 MODIS 两 种 植被 指数 相 比 最 明显 的 区 别 。 
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图 5 3 种 指数 拟 合 前 后 残 差 箱 线 图 ( 只 保留 生长 季 (4 一 9 A) 高 质量 数据 点 ) 


Fig.5 Residual box chart of three indices after function fitting (only retain the high quality data points in the growing season(4—9 months) ) 


NDVI/EVI 


SIF/(mW m sr! nm! ) 


2007 2008 2009 2010 2011 2012 2013 
年 份 Year 


图 6 2007—2013 年 3 种 指数 时 间 序 列 曲线 对 比 


Fig.6 ”Comparison of three time series curve for 2007 一 2013 


计算 3 种 指数 2007—2013 年 间 7 年 时 间 序 列 数据 均值 得 到 图 7 中 年 均值 时 间 序 列 曲线 ,NDVI 明显 区 别 
于 另外 两 种 数据 ,在 3 月 底 便 进 入 快速 增长 阶段 ,直至 8 月 初 才 逐 步 开 始 衰减 ,形成 平坦 且 较 长 的 生长 季 峰 期 
形态 ,10 月 初 逐 渐 下 降 并 趋 于 平稳 。 而 EV 与 SIF 曲线 的 变化 趋势 基本 一 致 , 均 在 4 月 未 开始 快速 上 升 ,但 
SIF 曲线 在 6 月 中 旬 便 较 早 到 达 生 长 季 峰 值 ,之 后 迅速 衰减 ,两 者 曲线 在 7 月 中 旬 逐 渐 重 合 ,EVI 稍 滞后 于 
STF ,9 月 未 结束 衰减 过 程 并 趋 于 稳定 。 
3.2 3 种 指数 物候 特征 参数 对 比分 析 

3 种 数据 在 时 间 序 列 曲线 形态 上 有 较 大 差别 ,得 到 的 物候 参数 差异 较 大 。 如 表 2 所 示 , 生 长 季 的 开始 时 
间 上 ,NDVI 早 于 另 两 种 数据 ,在 部 分 年 份 比 SIF 早 近 1 个 月 的 时 间 ,EVI 与 NDVI 相差 较 小 , 只 在 2007 年 有 近 
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16 d 的 差别 ,其 余年 份 基本 一 致 。 生 长 季 结 束 时 间 上 ， ore 
NDVI 平均 晚 于 SIF 1 个 月 左右 ,与 EVI 差 值 最 小 值 在 6 
d 左右 ,最 大 值 (2013 年 ) 达 到 了 25 d, 在 生长 季 的 持 
续 时 间 上 ,同样 是 NDVI>EVI>SIF ,在 部 分 年 份 NDVI 得 
到 的 LOS 近 2 F SIF 得 到 的 LOS, 显 示 出 了 巨大 的 差 
异性 。 

由 图 8 可 以 清晰 看 出 ,SIF 数据 的 季节 峰值 时 间 点 
明显 早 于 MODIS 数据 ,同时 也 早 于 SIF 数据 的 季节 中 
点 ,而 MODIS 数据 的 峰值 点 与 季节 中 点 则 更 为 接近 。 [2.3 456 7 8 9 10 11 12” 

由 于 季节 峰值 时 间 点 是 根据 每 一 期 数据 统计 得 到 ,因此 ay Moms 

受 限 于 数据 采集 的 时 间 跨 度 无 法 同 季 节 中 点 (计算 得 图 7 3 种 指数 7 年 均值 时 间 岸 列 旧 线 

出 ) 一 样 精确 ,因此 此 处 分 析 认 为 季节 峰值 时 间 点 与 季 Fig.7 The time series curye-mean values of three indices 
节 中 点 相差 在 一 个 数据 期 数 内 就 代表 时 间 基 本 一 致 。 averaged by seven years 

基于 这 种 判断 ,MODIS 数据 的 季节 中 点 与 季节 峰值 时 

间 点 基本 在 同一 时 间 , 而 SIF 数据 表现 为 季节 峰值 时 间 点 早 于 季节 中 点 约 2 一 3 个 数据 周期 。 


NDVI/ EVI 
SIF/(mW m? sr! nm”! ) 


R2 3 种 指数 的 物候 特征 参数 


Table 2 Phenological characteristics-of three indices 


年 份 NDVI EVI SIF 
Year SOS EOS LOS SOS EOS LOS SOS EOS LOS 
2007 98 287 189 115 278 163 132 270 138 
2008 108 288 180 115 274 159 122 275 153 
2009 117 280 163 122 271 149 116 269 153 
2010 122 295 173 125 285 160 127 278 151 
2011 121 287 166 127 273 146 143 263 120 
2012 113 300 187 113 287 174 126 269 143 
2013 108 302 194 110 276 166 128 259 131 


NDVI: 归 一 化 植被 指数 ,Normalized Difference Vegetation Index;EVI: 增 强 型 植被 指数 , Enhanced Vegetation Index; SIF; It 2k ¥é 56 , chlorophyll 
fluorescence ; SOS :生长 季 开 始 时 间 ,start of growth season; EOS :生长 季 结 束 时 间 ,end of growth season; LOS :生长 季 长 度 ,length of growth season 


以 上 提取 的 物候 参数 反映 出 了 NDVI、EVI 与 SIF 数据 物候 监测 结果 的 较 大 差异 ,植被 绿 度 在 叶绿素 还 未 
进入 活跃 期 时 便 已 开始 发 生 明 显 变 化 , 且 变 化 幅度 为 对 称 形态 。 而 叶绿素 的 高 强度 活动 周期 更 加 集中 和 短 
促 , 在 植被 进入 到 生长 季 之 后 迅速 达到 最 旺 戌 的 状态 ,而 后 快速 衰减 。 植 被 叶片 的 生长 也 在 叶绿素 光合 作用 
达到 峰值 后 进入 到 生长 最 快速 的 阶段 ,之 后 在 植被 仍 维持 一 定 的 绿 度 时 光合 作用 的 活跃 期 结 
3.3 3 种 指数 监测 森林 物候 能 力 的 比较 

利用 .2003 一 2010 年 间 7 年 间 通 量 站 点 的 CPP 逐日 数据 对 同时 段 内 的 NDVIEVI\SIF 得 到 的 物候 监测 结 
果 进 行 验证 。 

图 9 中 清晰 地 显示 出 3 种 指数 与 CPP 之 间 的 变化 关系 。NDVI 曲线 与 CPP 数据 差异 较 大 ,着 重 体 现在 
NDVI 监测 到 的 生长 季 明 显要 长 于 CPP 监测 的 生长 季 , 在 非 生长 季 与 CPP 数据 也 出 现 了 较 大 偏差 。EVI 相 较 
于 NDVI 则 有 较 大 改善 ,整体 上 较为 符合 GPP 数据 的 变化 趋势 ,但 非 生 长 季 数 值 明显 偏 高 ,无 法 捕捉 生长 季 
中 植被 快速 生长 时 的 变化 细节 。SIF 曲线 形态 与 CPP 之 间 存 在 良好 的 相关 性 ,不 仅 与 CPP 数据 点 基本 重合 ， 
上 且 在 生长 季 期 间 很 好 地 反映 出 植被 CPP 数值 的 变化 细节 ,总 体 偏 差 较 小 。 

图 10 中 ,4 种 数据 的 年 均值 曲线 图 也 反映 出 了 上 述 特点 ,SIF 曲线 依然 是 与 CPP 最 为 吻合 的 曲线 ,而 
MODIS 植被 指数 虽然 在 整体 趋势 上 也 保持 了 与 CPP 相同 的 变化 过 程 ,但 偏差 较 大 ,存在 着 约 1 个 月 左右 的 时 
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间 偏 差 。 综 合 图 9 图 10 中 的 曲线 对 比 关系 ,SIF 显示 出 PE 
与 植被 CPP 数据 最 好 的 匹配 关系 ,但 SIF 曲线 依然 出 一 SIF 峰 值 点 ”一 EVLA — NDVI 
MTEF GPP 开始 衰减 的 趋势 ,这 也 是 SIF 曲线 与 其 mis e 
他 数据 相 比 的 一 个 特性 。 2012 
2011 
mere : 
4 讨论 oe 
ks 
4.1 NDVI 与 EVI 时间 序列 曲线 的 差异 È 2009 


EVI 在 计算 公式 中 引用 蓝 波段 并 针对 残留 气 洲 胶 2008 
的 后 处 理 改 善 了 NDVI 的 不 足 ,使 两 种 指数 在 物候 曲线 2007 


产生 了 较 大 差异 ' 引 。NDVI 生成 算法 中 只 采用 了 红 光 160 170 180 190 200. 210 220 一 230 
和 近 红 外 两 种 波段 , 当 植 被 覆盖 度 较 高 时 ,植被 叶片 中 时 序 Day of year 
的 叶绿素 使 红 光 波段 很 快 趋 于 饱和 ,加 之 算法 本 号 为 非 图 8 3 种 指数 生长 季 峰 值 与 时 间 中 点 对 比 图 


线性 增长 浮 数 ,导致 高 植被 覆盖 度 时 出 现 饱 和 效应 , 低 Fig.8 The peak value in growing season and the midpoint of 
植被 覆盖 度 时 数值 易 偏 高 夸大 曲线 振幅 !55 ,因此 曲线 ree Indices 

呈现 出 较 早 开始 上 升 .生长 季 较 长 昌平 坦 ， 而 NDVI 后 图 中 曲线 为 3 种 数据 的 季节 峰值 时 间 点 的 拟 合 曲线 

处 理 中 对 大 气 校正 的 不 彻底 导致 水 汽 、 气 溶胶 残留 ， 


EVI SIF 
35 1.0 6 


NDVI/EVI 
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图 9 2007 一 2010 年 MODIS 植被 指数 (NDVI、EVI) SIF 时 间 序 列 曲线 与 GPP 对 比 图 
Fig.9 Comparison of MODIS vegetation index( NDVI, EVI) and SIF time series curve with GPP for 2007 一 2010 


此 在 长 自 山 站 点 春 冬 季 存 在 冰雪 和 履 盖 的 情况 下 ,更易 受到 云雾 迟 掩 ,冰雪 覆盖 影响 ,使 得 非 生长 季 期 间 出 
现 较 大 的 起 伏 以 及 多 处 异常 峰值 。 
4.2 SIF 与 植被 指数 的 差异 
相 比 反映 植被 绿 度 的 植被 指数 ,通过 探测 植被 光合 作用 发 射 的 严 光 信号 ,快速 .无 损伤 反映 植被 光合 作用 
特征 的 叶绿素 荧光 , 则 成 为 了 监测 植被 对 光 能 利用 情况 的 “ 探 针 ”同时 也 可 以 很 好 反映 植被 的 生长 状态 。 
SIF 与 植被 指数 之 间 的 差异 主要 体现 在 两 个 方面 。 在 本 研究 中 ,SIF 物候 曲线 明显 晚 于 两 种 植被 指数 进 
入 生长 季 , 且 较 早 结束 ,大 约 在 4 月 末 开 始 ,9 月 末 结 束 。 刘 新 杰 等 .对 我 国 西北 、 西 南 、 华 北 地 区 进行 物候 人 研 
究 的 结果 显示 SIF 高 值 通常 在 4 一 5 月 期 间 , 而 MODIS 数据 高 值 通常 出 现在 7 一 9 月 。 本 文中 SIF 的 高 值 出 现 
时 间 大 致 为 6 一 7 月 ,植被 指数 高 值 出 现时 间 大 致 为 7 一 9 月 ,与 刘 新 杰 等 的 研究 具有 相似 的 规律 性 , 即 SIF 要 
FPF MODIS 植被 指数 达到 生长 季 峰 期 水 平 ,但 SIF 峰 期 提前 的 时 间 与 本 文 存在 1 一 2 个 月 的 差异 。 考 虑 到 该 
研究 中 的 3 处 区 域 中 ,西北 区 与 本 文 差 值 最 小 ,西南 最 大 ,因此 水 热 条 件 越 充 足 , 植 被 的 SIF 曲线 到 达 峰 期 的 
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图 10 4 种 数据 年 均值 时 间 序列 曲线 图 
Fig.10 Time series curves of mean value of NDVI, EVI, SIF and GPP 


时 间 就 越 早 。Walther 等 ”对 北美 中 高 纬 落叶 林 进 行 物 候 研 究 并 对 比 了 MODIS NDVI, EVI, GOME-2 SIF 数 
据 , 研 究 结果 与 本 文 总 体 一 致 ,SIF 在 4 月初 快速 增长 ,5 月 末 达 到 最 大 值 ,7 一 9 月 缓慢 下 降 并 逐步 进入 秋季 
快速 衰减 阶段 。 植 被 指数 的 EOS 晚 于 SIF 的 情况 ,通常 是 反映 植被 绿 度 变 化 的 指数 共有 的 问题 ,因为 在 植被 
进入 秋季 衰老 期 后 ,植被 光合 作用 虽然 大 幅度 减弱 并 趋 于 停 i 但 叶片 绿 度 并 不 会 迅速 反映 这 种 改变 而 是 存 
在 一 个 渐变 过 程 ,日 卫星 信号 通常 会 被 衰老 死亡 但 未 凋落 的 叶片 影响 ,因此 相 较 于 SIF 的 秋季 衰减 NDVIEVI 
均 存 在 一 个 滞后 过 程 汪 。 

SIF 在 夏季 生长 季 会 出 现 早 于 植被 指数 的 快速 下 降 , Yang 等 :中 利用 摄像 机 与 MODIS 图 像 进行 林 冠 测量 
时 提 到 过 这 种 现象 , Walther 等 利用 GOME-2 反 演 的 SIF 数据 对 中 高 纬 森林 进行 物候 研究 时 也 发 现 北美 的 
多 处 落叶 阔 叶 林 以 及 部 分 混交 林 出 现 这 种 情况 。 产 生 这 种 现象 的 原因 ,Wilson 等 “提出 夏季 和气 元 素 分 配 的 
改变 导致 光合 作用 中 一 种 可 能 有 重要 作用 的 Rubisco 酶 发 生 了 改变 。Bauerle 等 ”通过 实验 提出 是 光 周 期 的 
调节 使 得 光合 作用 在 夏季 又 降 。 本 文 研 究 中 发 现 SIF 曲线 夏季 又 降 只 在 落叶 林 以 及 混交 林 中 出 现 , 且 时 间 点 
在 夏至 日 前 后 对 于 北半球 的 中 高 纬 地 区 来 说 ,夏至 日 前 后 太阳 高 度 角 达 到 一 年 中 最 大 值 ,气温 开始 升 
高 ,气候 学 上 以 每 5d 的 来 均 气温 高 于 22°C 的 始 日 作为 夏季 开始 ,并 伴 有 频繁 降雨 天 气 “ 。 综 合 以 上 两 点 ， 
落叶 林 ( 以 及 如 本 文中 具有 明显 落叶 林 特 性 的 混交 林 ) 的 生长 习性 使 其 在 春季 复苏 期 会 吸收 更 多 的 二 氧化 碳 
进行 光合 作用 以 生出 大 量 新 叶 “| ,在 这 过 程 中 光合 作用 强度 迅速 增 大 。 但 随 着 夏至 日 前 后 光照 强度 与 时 长 
的 增加 ~ 温度 的 升 高 ,使 叶绿素 的 光合 作用 受到 抑制 “ ,因此 伴随 光合 作用 强度 的 下 降 SIF 也 表现 出 在 夏季 
迅速 衰减 的 特性 。 
4.3 3 种 指数 与 通 量 站 GPP 数据 的 关系 

GPP 是 表征 植被 通过 光合 作用 产生 有 机 物 与 能 量 总 和 的 生物 变量 , 通 量 塔 高 时 间 频 率 的 日 观测 数据 计 
算得 到 的 GPP 数据 对 植被 的 生长 状态 具有 非常 高 的 代表 性 。GPP 与 SIF 之 间 的 季节 变化 相关 性 已 在 多 
种 植被 类 型 上 被 验证 449 | Yang 等 (更 是 指出 SIF 数据 ,包括 卫星 与 地 面 两 种 获取 方式 ,在 落叶 林 类 型 上 
与 通 量 塔 CPP 存在 每 日 和 季节 性 的 强 关 联 。 本 文 的 研究 中 SIF 曲线 虽然 存在 夏季 又 降 现 象 ,但 依然 显示 出 
了 与 通 量 站 GPP 数据 的 良好 一 致 性 ,Bauerle" 等 通过 试验 证 实 这 种 变化 与 叶片 光合 作用 的 季节 变化 一 致 。 
这 种 下 降 并 未 使 得 CPP 出 现 明 显 的 下 降 , 说 明 虽 然 光 合作 用 的 强度 受到 抑制 ,但 产生 的 有 机 物 总 量 却 没有 受 
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到 较 大 干扰 ,考虑 到 SIP 强度 与 光 能 利用 率 (LUE) 之 间 存 在 可 靠 负 相关 关系 ,因此 本 文 推测 SIF 虽然 快速 
下 降 , 但 叶绿素 利用 太阳 辐射 进行 光合 作用 的 效率 却 得 到 了 提升 ,因此 GPP 并 未 出 现 明显 下 降 。 植 被 指数 的 
局 限 性 导致 与 CPP 之 间 存 在 明显 偏差 ,EVI 相 比 于 NDVI 与 GPP 吻合 程度 更 高 一 些 ,虽然 无 法 很 好 捕捉 植被 
生长 变化 的 特征 细节 ,但 对 于 季 闻 性 变化 较 大 的 和 森林 类 型 ,凭借 叶片 的 物候 特征 与 二 氧化 碳 吸收 具 有 强 关 联 ， 
因此 也 可 以 用 来 监测 GPP AER 。 


5 结论 


本 文 以 长 白山 针 阔 混交 林 站 点 为 研究 区 域 ,采用 MODIS NDVI EVI 与 GOME-2 SIF 数据 ,利用 数据 拟 合 
与 动态 阐 值 法 提取 物候 特征 参数 ,并 对 比分 析 3 种 数据 在 森林 植被 物候 监测 上 的 区 别 与 各 自 特性 ,最 后 利用 
通 量 站 GPP 数据 进行 验证 。 

研究 结果 表明 ,NDVI 受 限 于 计算 公式 和 后 处 理 中 的 不 足 使 得 噪声 影响 明显 ,在 非 生 长 季 数 据 起 伏 波动 
最 大 ,生长 季 出 现 饱和 现象 ,与 CPP 数据 也 存在 较 大 差异 。EVI 曲线 更 加 平滑 规律 且 与 CPP 二 配 程度 更 好 ， 
但 秋季 衰减 时 间 稍 晚 于 GPP 数据 。SIF 虽然 存在 波动 但 总 体 较 小 ,季节 变化 趋势 明显 ,与 CPP 季节 变化 基本 
一 致 , 较 好 反映 了 植被 生长 的 变化 特征 ,但 受 夏 至 日 前 后 光照 强度 和 温度 升 高 影响 出 现 夏 季 又 降 现 象 。 

结合 以 上 人 研究 结果 ,植被 指数 对 于 植被 绿 度 的 监测 虽然 凭借 生长 过 程 中 与 三 氧化 碳 通 量变 化 的 关联 可 以 
一 定 程度 上 监测 GPP 的 季节 变化 反映 植被 的 生长 状态 ,但 受 限于 计算 与 后 处 理 方式 .植被 的 季节 变化 幅度 
等 ,在 关键 物候 时 间 上 普遍 存在 提前 和 光 后 现象 。SIF 与 CPP 上 共有 良好 二 致 性 ,可 以 较 好 反映 植被 的 生长 过 
程 并 精确 捕捉 植被 光合 作用 的 变化 细节 。 但 反 演 SIF 菊 光 波段 的 选择 ,后 处 理 中 质量 控制 的 方法 都 会 影响 数 
据 在 物候 研究 中 的 参考 性 ,并且 SIF 的 变化 与 多 种 植被 参量 (APAR ,LUE 等 ) 的 关系 依然 不 够 明确 ,因此 对 于 
SIF 变化 细节 的 研究 依然 是 吸 待 解决 的 问题 。 
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